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I ■ Abstract 

The effects of spatially correlated noise on a phenomenological 
equation equivalent to a non-local version of the Kardar-Parisi-Zhang 



Q ' (KPZ) equation are studied via the dynamic renormalization group 

^ ■ (DRG) techniques. The correlated noise coupled with the long ranged 

nature of interactions prove the existence of different phases in differ- 

f^ ent regimes, giving rise to a range of roughness exponents defined by 

t:;;^- ■ their corresponding critical dimensions. Finally self-consistent mode 

On ! analysis is employed to compare the non-KPZ exponents obtained as a 

_i I result of the long range-long range interactions with the DRG results. 
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The interests in the non-equihbrium growth mechanisms in the forma- 
tion of surfaces and interfaces as also the description of directed polymers, 
bacterial growth, etc. and of late the protein folding problems, although 
apparently being the representation of different physical processes, have all 
been encapsulated in one single non-linear continuum equation, the much cel- 
ebrated Kardar-Parisi Zhang equation [1] . The notion of universality classes 
defined by this standard KPZ equation [1] suggests the existence of a phase 
transition from the EW class [2] to the non-linear KPZ class above a particu- 
lar critical dimension {d > 2). However although this stochastic equation has 
by now become a model of dynamic critical phenomena for a vast range of 
growth problems, the fact remains that the basic non-linearity studied here 
is of a local nature, that is to say the growth occurs along a continuously 
varying local normal. 

To incorporate the long-ranged nature of interactions which are necessary 
for a wide class of problems, eg. the long-ranged hydrodynamic interactions 
[12-14], proteins, colloids [15-18], etc. SM and SMB [3] developed a Langevin 
type equation studying the effects of the long-ranged (LR) feature of an evolv- 
ing surface going beyond the local description of the KPZ non-linearity for 
the case of white noise. There the approach essentially consisted in introduc- 
ing a term in the basic Langevin equation capable of correlating each site of 
the growing surface with all other sites. The objective was the transforma- 
tion of the local non-linear term representing the lateral growth beyond the 
strict local description such that the correlation length now becomes at least 
of the system size. Still there remains to be seen the effects of interaction of 
correlated colored noise with the KPZ or KPZ-type non-linearity. 

The results of non-white noise for the growth of rough interface has been 
generalized by Medina et al [4]. In two remarkable papers Chekhlov and 
Yakhot [5,6] and Hayot and Jayaprakash [7] have observed the effects of 
correlated noise for the one-dimensional Burgers' equation. They explored 
the occurence of shocks as well as the large distance, long time statistics of 
the fluctuations. Working along this line Frey, Tauber and Janssen [8] have 
also resorted to the long-ranged description of noise in their treatment of the 
scaling regimes and critical dimensions in the standard KPZ problem as well 
as in the conserved case [9]. Actually in all these cases, the nature of the 
noise is determined from the fact that to maintain turbulence in the flow, 
energy has to be supplied at large length scales near the boundaries and the 
consequent Kolmogorov type of dimensional argument brings about a spatial 
dependence in the noise correlation. 

Starting with the non-local equation proposed in [3], we have gone one 
step further in putting forth the effects of a non-white spatially correlated 
noise akin to that used in [4], the objective being the analysis of special 
features of interaction of long-ranged nature of KPZ type non-linearity with 
the long-ranged spatial correlation in noise. Use is made of dynamic renor- 
malization group (DRG) techniques in arriving at the dynamic exponents, 
etc. We see that even for d < dc, both weak and strong noise interacting 



both with the local and non-local nature of the non-linearity give a range of 
critical exponents spanning a four dimensional space in terms of the dimen- 
sionless interaction strengths. Finally we reconfirm our DRG results from 
the self-consistent mode power counting arguments in the line developed in 
[10,11]. 

The starting point of our analysis is the equation: 



dh(f t) 1 

^ ' ^ - z/V^/i(f,t) + - / drv{r')Vh{f+r',t).Vh{f-r',t)+r]{f,t) (1) 



dt ' ' ' ' 2 

where z/ is the diffusion constant and i]{f,t) is the noise defined by 

< r]{k,uj)r]{k',uj') >= 2D{k)5'^(k + k')5{uj + uj') (2) 

Again going by the prescription of [3], the kernel v{r) has a short-ranged 
part Ao(5'^(r) and a long-ranged part r''"'^. In Fourier space, 

v{k) = Ao + \pk-P (3) 

All standard KPZ results are expected for \p = 0. Our introduction of the 
long-ranged (LR) version of the noise in coupling with the long-ranged (LR) 
part in the interaction spectrum is expected to define a legion of continuously 
varying exponents for a set of new non-KPZ fixed points. 

The two exponents of interest, the roughness exponent a and the dynamic 
exponent z come along with the two-point height correlation function in the 
hydrodynamic limit {k,uj) -^ 0, 



< h{k,uj)h{k',uj') >^ \kf '^ ^"^ 6'^{k + k')6{uj + uj')f{^^) (4) 

I k I 

All information regarding the dynamic universality class of the phase will 
be contained in this a and z. 

At d=l, the values of a{= |) and z{= |) can be exactly determined. But 
at d=2, there is a transition from the Gaussian fixed point (EW) and the non- 
linearity grows under rescaling. Simple scaling from x -^ bx, h ^> h°'h and 
t -^ b^t shows that both the short (SR) and long-ranged (LR) contributions 
in the interaction kernel are relevant for rf < 2 (where by SR interaction we 
mean the standard KPZ type non-linearity, that is the kernel is a constant 
multiplying 6'^{r) and the LR interaction implies a non-KPZ type r dependent 
part) . Under this scale transformation the parameters of the equation change 
by z/ ^ 6^-V, Ao -^ b'^+'-^\o,\p -^ b'^+'+f'-'^Xp. If the noise strength D{k) 
in the hydrodynamic limit {k, a; ^ 0) is given by D{k) = Dq + D^rk^'^'^, then 
Do -^ b-'^'^'^+'Do and D„ -^ 6-4a-d-2a+^^^^ For 2 < rf < 4 = 2+2p-4a, Ap 
is relevant at the EW fixed point (z=2) and for p > 0, a non-KPZ fixed 
spectrum should be the outcome. The following DRG analysis gives a horizon 
of unfounded results which shrink to the known results in d=l and 2 as in 



[3], new stable points are obtained and furthermore the introduction of the 

non-local noise provides even more complexity in the interaction. 

^ ah(k,Lu) ^ 

The generalized Green's function G(k, u) = ,.,^^''i'°''"''' — 7- in the sym- 
metrized Fourier space is given by, 

1 r rl'^k' duj' 

G{k,u;) = GoiKoo) + 4.{-fGo\k,u;) J -^ — k'^.k'_ 

v{2k') v{— + k')G,{k'_,uj'_) Goik'^, co') 
Go{-k'^,-uj'^)2Dik'^,Lu'^) (5) 

where symbolically X± = y ±X', with X = poi u. Go^k, u) = —iu + uk"^ 
represents the unperturbed term in the propagator expansion. We consider 
diagrams upto the one- loop order 0{v'^) and follow this expansion with a 
detailed dynamic renormalization procedure as in Medina, et al [4]. The 
propagator expansion gives us the viscosity flow, the vertex renormalization 
gives the flows of Aq and Xp and the noise renormalization gives the Dq and 
Dcr flows. As usual we integrate out the fast modes in the momentum shell 
Ae~' < /c < A, set the cut-off A = 1 and develop the flow equations of the 
rescaled parameters around the fixed points. 

Due to the Galilean invariance of eqn.(l), Aq is not renormalized. The 
RG transformation being analytic in nature, Xp is also not renormalized, only 
the Galilean identity is modified in this case (2-p instead of 2), 

a + z = 2-p (6) 

where p = for Aq fiow. From the above considerations, we get the 
following fiow equations for z/, D's and X's, 

- = u[z-2-K, ] (7) 

^ = Dik) [z-2a-d- m] +K,^ v\2) (8) 

-^ = X^[a + z-2 + x\ (9) 

where x = or p respectively. Here /(g) = ^^^^ ' \k=q and g{q) = 

Q^^^' \k=q) Kd = -i0rd., Sd representing the d-dimensional surface of a unit 
d+1 dimensional sphere. 

In terms of the dimensionless interaction strengths Uq^ = Kd ^° f'" , where 
short- ranged interaction couples with long-ranged noise (the Medina, et al 
zone) and similar other parameters U"^^^ = Kd ^^3 "^ , Uqq = Kd '^°^f " (ordi- 
nary KPZ case), UpQ = Kd ''^3 ° (the SM and SMB zone), the fiow equations 
can be combined to give. 



dUo,a 2-d + 2a d-2-2a 

-IF - f 2 J^o- + 3i ]f/o,. + 

^ [3(rf - 2)Ul, + 3.2-^(rf - 2 - 3p)f/,^,o + 

3(1 + 2-0 (rf - 2)f/o,of/p,o + 3.2-^(rf - 2 - 2a - 3p)f/2^ 

+ 3(l + 2-0(c?-2-2a)[/o,.f/p,a] (10) 



rff/p,. _ 2-rf + 2a + 2p ^_^^ rf-2-2a-3p ,^^3 



^ [3(rf - 2)f/Jo + 3.2-''(rf - 2 - 3p)f/,% + 

3(l + 2-'')(rf-2)f/o,of/p,o + 3.(rf-2-2(T)?7o'_, + 

?,{l + 2-P){d-2-2a)Uo,.Up^„] (11) 






^[2-''{(3 + 2-'')rf-6-9p}f/Jo 

+ {3(1 + 2-P){d - 2) + rf.2-''+^}f/o,of/p,o + 3(rf - 2 - 2a)f/2, 
+ 3.2-^(ci - 2 - 2cr - 3p)f/2^ + 3(1 + 2-f>){d - 2 - 2a)f/o,.f(A2) 



^r - ^ — 2 — J^^'° + ^ M J^^'° 



^[(4rf-6)f/o% + {3(l + 2-0(rf-2) + rf.2-^+i}f/o,of/p,o + 

3(rf - 2 - 2o-)Ul^ + 3.2-^(rf - 2 - 2ct - 3p)f/^2^^ + 

3(l + 2"0(rf-2-2a)f/o,.f/p,.] (13) 

Now let us define two sets of parameters: Rq = 77^, R^ = jr^ and Sq = 

Uo,o q Up,o 

Uo.a ' '^P ~ Up,„ ■ 

We see that -^ = —pRy, wliere y=0, a and consequently this rules out 
any off-axial fixed point in the Rq, R^, parameter space (except for the trivial 
(7 = case). 

In the ?7o,o, f^p,o plane, the axial fixed points are given by, 

SR-SR^ U*\,o = 0, U*\,o = ^^^^ a + z = 2 (14) 



ro CD— rr*2 n rr*2 ^v*^ ^ 2p) 

^^-^^=^ 0,0 = 0, f/ P,o= ^.,^(3^^_,^^_g_g^^ ,c. + . = 2-p. 



(15) 



The first set (SR-SR) gives the well-known KPZ fixed point with z = 
3/2, a = 1/2 for d=l. But the second set (LR-LR)gives the non-KPZ be- 
haviour and the results exactly match with [3] in this Uqa, Up^ plane: 



_ (rf-2-2p)(rf-2-3p) 

^lc/*\o=o - ^+ p + 2-P)-6-9p] ' 

« 1^-0,0=0 = -P- [(3 + 2-P) - 6 - 9p] ^^^^ 

In diagram B of Fig.l, the dotted line gives an unstable zone between 
d = TT^ to d = 2 + 2p. Above the critical dimension, diagram C shows a 
smooth phase. For p = 0, all the LR fixed points go over to the SR ones. 

In the ^0,0- , Up^o- plane defined by R^, there also are only two sets of axial 
fixed points: 

SR-LR= U*\,^ = 0, U*\,^ = y (17) 

.<? .9 Ad(2 - d + 2a + 2p) , , 



3.2~p{2 + 2(r + 3p-d)' ^ ' 

In this case the phase diagrams are exactly identical to Fig.l, only the 
critical dimension is now modified to rfc = 2 + 2p — 4o" and the unstable zone 
now lies between c? = 2 + 2o" + 2p to c? = 2 + 2(7 + 3p in this plane: 



^lc/*V=o = 3(4 + rf-2a-2p), 

a|c/'V=o = -p+^(2-rf + 2a + 2p). (19) 

However a completely different qualitative behaviour is observed with 
Sq and S^. The Sq fiow equation reads ^ = Sq{-g + \UIq + 2~2pf/2^ + 

^t/o,of^p,o and as such off- axial fixed points exist in this case for the 4- 
dimensional space of U's. But in the f/o^O; ^o.o- plane we get only axial fixed 
points where now f/o,o? f^p,o =constants. Here the axial fixed points are 

(20) 

(21) 
This gives 
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(22) 



The results exactly match with the Medina, et al predictions and the 
phase diagram is given by FIG. 2, X and Y are the axial and off-axial fixed 
points respectively. The point to be noted here is that unlike the previous 
two cases there is no unstable excluded region in the non-KPZ case. 

Although similar arguments as above apply for the Sp flow, the most 
spectacular results are seen in the f/p,o, Up^a plane where the fixed points are 
given by. 



.2 ^ . rr.2 ^ 4d{2-d+2p) 

"'" ' ^'° 2-p[Q + 9p - (3 + 2-p)] 



LR-SR^ U*\,^ = 0, f/*%,o = .,.: : I'' (23) 



.2 ^ n U*^ = ^d{2-d + 2a + 2p) 
"'^ ' "'" 3.2-^(2 + 2a + 3p-rf)' 



LR-LR^ w\, = 0, w\,^ = ^7::: : : '\ - m 



Here for both the sets we have unstable regions bounded by 2 + 2p > d > 
1^ (LR-SR) and 2 + 2a + 2p < d < 2 + 2a + 3p for p, a > (LR-LR). 
Also both these fixed points give non-KPZ results: 



1 
3 
and 



z\u*%,o=o= ;j(c^ + 4-2a-2p) (25) 



. I - 9 M id-2-2p)id-2-3p) 

z \u*^,^=o - 2 + (3 ^ 2-P)ci - 6 - 9p ^^^^ 

and as such the phase diagram in Aq, Xp, Dq, D„ is actually on a 4- 
dimensional space with both axial and non-axial fixed points. 

Now use is made of self-consistent mode analysis (exact in the spherical 
limit) to generate the non-KPZ exponents in this complex space where both 
the SR-LR noises are interacting with the SR-LR non-linearity. 

The Dyson's equation gives G^^{k,uj) = —iuj + uk"^ + S(A;,c^), where 
S(A;,a') is the self-energy term and the corresponding scheme gives 



S(^,a;) = / '^'^k'^.k'^v{2k')v(^ + k') 



d'^k' duj' p p p 3k 

— J- k:.k v{2k )v{ — 

271'^ 2n + ~ ^ ' ^ 2 

Go{k'_,Lj'_) Go{k'_^,uj') Go{-k'^, -uj'_^) 2D{k'^,u'^) (27) 
The noise correlation function at the same order of accuracy (0(f ^)) gives 

d^k' du' 



D{k,uj) = Do{k,uj)+ \Go\k,uj)\ J 



(27r)°' 2tt 

k'^\k'J v\2k') D{k'^, uj'^) D{k'_,uji) 
Go{k'^,uj'_)Go{-k'_,-u'_)Go{k'_^,u;'_^)Go{-k'^,-u;'_^) (2J 

We now make the following scaling ansatz for S(A;,a;) and D{k,uj): 



J:{k,uj) = \kfij{-^) (29) 

I ^ I 

D{k,u;) = \kf'''<l){-^) (30) 

I ^ I 
In eqn.(20) the noise strength Dq will be renormalized only if the second 
term is more singular than the first. If the second term in eqn.(20) dominates 
in the long wavelength limit, eqn.(20) reduces to 

v\2k') D{k'^,uj'^) D{k'_,ujL) 

G{k'_,uj'_) G{-k'_, -u'_) G(k'^, uj'^) G{-k\, -uj'^) (31) 

In the limit a; — * 0, taking D[k, 0) = D„ k^^" and v{k) = Xp k~'', simple 
power counting from eqns.(19) and (23) gives 

1 1 

z = -{d + A-2a-2p), a = -p + - {2 - d + 2a + 2p). (32) 
o o 

This values of z and a are seen to tally very well with our dynamic RG 
derivations in the Uofi, f/p,o, f^o,o-, f^p.o- space, where in either of the two planes 
the non-KPZ values coincide with eqn.(24). However, the Upfi, Up^a plane 
is special in that it provides two axial non-KPZ values only one of which 
appears in the self-consistent results. The other value is actually an off- 
shoot of the SM and SMB derivation [3] where the relevant LR non-linearity 
is interacting with the SR part of the noise. That way we are trodding into a 
new domain of roughness where both the non-KPZ values of z (and a) exist 
although no rough-rough phase transition apparently takes place. 

In summary, we have started with a simple phenomenological equation 
which incorporates a non-local term in its interaction spectrum and while 
coupling with spatially correlated noise develops a set of dynamic and growth 
exponents which contain a non-KPZ part. As and when the SR part in the 
spectrum becomes unstable, these non-KPZ domains surface up and com- 
pletely different critical behaviour sets into existence. With negative values 
for the long and short-ranged non-linearities the phase diagrams are mod- 
ified, with the LR roughness now giving way to SR roughness without the 
appearance of any excluded instability in the phases, the only exception be- 
ing in the case of the special plane already discussed. Also it would be 
worthwhile to mention that although the non-local contribution in the non- 
linearity never generates it's short-ranged counterpart, the non-local part 
in the noise spectrum develops a white noise. We reconfirm all these DRG 
observations from a self-consistent technique and arrive at the same set of 
non-KPZ exponents when the noise strength remains non-renormalized. Fi- 
nally comparisons with established results [3] and [4], which constitute only 
parts of our whole domain, provide expected results. 
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2 Figure Captions 

Fig.l \p vs Aq phase diagram. The solid lines along x-axis represent LR 
phases while the dotted line in (B) shows a smooth phase. 
Fig. 2 Dq vs Dfj flow In the f/o,o vs f/o,o- plane. X and Y are the axial and 
off-axial flxed points respectively. 



10 



X 



I 



-p I (A) d<2 
LR rough 



Xp I (B) 2+2p-4|3>d>2 



Xr, 



LR rough 



■ An 



LR rough 

• • ^^^^ 



(C) d>2+2D-43 
LR rough 



■ Xn 



LR rough 



LR rough 



Figure L XpVS Xo phase diagram. The solid lines along the X-axis represent 

the LR rough phase while the dotted line in (B) represents a smooth 
phase for p > 0. 
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Figure 2. Do vs Da flow diagram. X is the axial fixed point and Y gives the off-axial 
fixed point. 



